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ABSTRACT 

The collapse and fragmentation of initially prolate and oblate, magnetic 
molecular clouds is calculated in three dimensions with a gravitational, radiative 
hydrodynamics code. The code includes magnetic field effects in an approximate 
manner: magnetic pressure, tension, braking, and ambipolar diffusion are all 
modelled. The parameters varied for both the initially prolate and oblate clouds 
are the initial degree of central concentration of the radial density profile, the 
initial angular velocity, and the efficiency of magnetic braking (represented by 
a factor fmb = 10~^ or 10"'^). The oblate cores all collapse to form rings that 
might be susceptible to fragmentation into multiple systems. The outcome of 
the collapse of the prolate cores depends strongly on the initial density profile. 
Prolate cores with central densities 20 times higher than their boundary densi- 
ties collapse and fragment into binary or quadruple systems, whereas cores with 
central densities 100 times higher collapse to form single protostars embedded in 
bars. The inclusion of magnetic braking is able to stifle protostellar fragmenta- 
tion in the latter set of models, as when identical models were calculated without 
magnetic braking (Boss 2002), those cores fragmented into binary protostars. 
These models demonstrate the importance of including magnetic fields in studies 
of protostellar collapse and fragmentation, and suggest that even when magnetic 
fields are included, fragmentation into binary and multiple systems remains as a 
possible outcome of protostellar collapse. 



Subject headings: hydrodynamics — ISM: clouds — ISM: kinematics and dy- 
namics — MHD — stars: formation 
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1. Introduction 

Fragmentation during protostellar collapse is widely accepted to be the primary mecha- 
nism for the formation of binary and multiple star systems (e.g., Lafreniere et al. 2008; Chen 
et al. 2008). While it is clear that the overall form of the initial mass function for stars is 
directly tied to the initial conditions for protostellar collapse, i.e., the mass function of dense 
cloud cores (e.g., Dib et al. 2008; Swift & WiUiams 2008), fragmentation is necessary for 
producing binary star systems within these individual dense cores (Lafreniere et al. 2008; 
Chen et al. 2008). 

Three dimensional calculations of the collapse of centrally condensed, rotating cloud 
cores have been computed for quite some time (e.g.. Boss 1993) and continue to attract 
theoretical attention (e.g., Saigo et al. 2008; Machida 2008; Commercon et al. 2008). 

These calculations neglected the effects of magnetic fields. However, observations of OH 
Zeeman splitting in dark cloud cores have shown that magnetic fields are often an important 
contributer to cloud support against collapse for densities in the range of 10'^ — 10"^ cm^'^ 
(Troland & Crutcher 2008). Given this observational constraint, it is clear that three dimen- 
sional hydrodynamical collapse calculations should include magnetic field effects as well as 
self-gravity and radiative transfer (e.g.. Boss 1997, 1999, 2002, 2005, 2007). Magnetic fields 
are now being included in other three dimensional collapse models as well (e.g., Machida et 
al. 2004, 2005a,b, 2007, 2008; Kudoh et al. 2007; Price & Bate 2007, 2008). In particular. 
Price & Bate (2007) found that while magnetic pressure acts to resist fragmentation dur- 
ing collapse, magnetic tension can actually aid fragmentation, confirming the results found 
by Boss (2002). Machida et al. (2004, 2005a,b, 2007, 2008) generally found that binary 
fragmentation could still occur provided that the initial magnetic cloud core rotated fast 
enough. 

Magnetic braking is effective at reducing cloud rotation rates during the pre-coUapse 
cloud phase, but has relatively little effect during the collapse phase, according to the two 
dimensional magnetohydrodynamics models of Basu & Mouschovias (1994, 1995a,b). How- 
ever, Hosking & Whitworth (2004) found that rotationally-driven fragmentation could be 
halted by magnetic braking during the collapse phase. Boss (2004) argued that the ther- 
modynamical treatment employed by Hosking & Whitworth (2004) could have been more 
important for stifling fragmentation than magnetic braking, but did not offer any models of 
magnetic braking to support this assertion. 

Price & Bate (2007) presented models of the collapse of magnetic cloud cores, finding 
that magnetic pressure was more important for inhibiting fragmentation than either magnetic 
tension or braking, contrary to the results presented by Hosking & Whitworth (2004) and 
Promang et al. (2006), who found no evidence at all for the fragmentation of magnetic clouds. 
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Promang et al. (2006) assumed ideal magnetohydrodynamics (MHD), i.e., without ambipolar 
diffusion, a fact that is hkely to stifle fragmentation, whereas Hosking & Whitworh (2004), 
Boss (2004), and Price & Bate (2007) included ambipolar diffusion. However, subsequent 
ideal MHD collapse calculations by Hennebelle & Fromang (2008) and Hennebelle & Teyssier 
(2008) found that magnetic clouds could fragment if the initial density perturbation was large 
enough, and they speculated on what would happen if ambipolar diffusion was included 
in their models. These fundamental differences in the results of magnetic cloud collapse 
calculations results highlight the need to compare models where only one parameter at a time 
is being changed, so that the true effect of changing that one parameter can be discerned. 
Such a comparison is the major goal of the present study. 

Banerjee & Pudritz (2006) considered the collapse of magnetized cloud cores, finding 
that even though considerable angular momentum was lost from the cloud core by magnetic 
outflows from the disk's surface, the cloud was still able to collapse and fragment into a close 
binary protostar system. 

Price & Bate (2007) and Hennebelle & Teyssier (2008) both considered the collapse of 
spherical, magnetic cloud cores with initially uniform density and uniform magnetic field 
strengths. Machida et al. (2004, 2005a,b) studied the collapse of initially cylindrical cloud 
cores in hydrostatic equilibrium. Machida et al. (2008) considered the collapse of cloud 
cores with density profiles appropriate for Bonnor-Ebert isothermal spheres, similar to the 
Gaussian radial density profile clouds studied by Boss (1997, 2002) and by this paper. Such 
centrally- condensed density profiles represent the best guesses as to the radial structure of 
pre-coUapse molecular cloud cores (e.g., Myers et al. 1991; Ward- Thompson, Motte, & Andre 
1999). 

Boss (2002) modelled the collapse of initially prolate and oblate cores, including several 
magnetic field effects, but ignoring magnetic braking. Prolate and oblate cloud shapes have 
been inferred from observations of suspected pre-coUapse molecular cloud cores (e.g., Jones, 
Basu, & Dubinski 2001; Curry & Stabler 2001). Here we use the magnetic braking ap- 
proximation developed by Boss (2007), originally applied to filamentary clouds, to examine 
the importance of magnetic braking for the same prolate and oblate models as those calcu- 
lated by Boss (2002), allowing a direct comparison between identical protostellar collapse 
models with and without magnetic braking. These models thus directly address the differ- 
ent outcomes of the models with magnetic braking but without detailed thermodynamics by 
Hosking & Whitworth (2004) compared to those without magnetic braking but with detailed 
thermodynamics of Boss (2004): the present models include both effects. 
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2. 



Numerical Methods 



The numerical models are calculated with a three-dimensional hydrodynamics code that 
calculates finite-difference solutions of the equations of radiative transfer, hydrodynamics, 
and gravitation for a compressible fluid (Boss & Myhill 1992). The hydrodynamic equations 
are solved in conservation law form on a contracting spherical coordinate grid, subject to 
constant volume boundary conditions on the spherical boundary. The code is second-order- 
accurate in both space and time, with the van Leer-type hydrodynamical fluxes having been 
modifled to improve stability (Boss 1997). Artiflcial viscosity is not employed. Radiative 
transfer is handled in the Eddington approximation, including detailed equations of state 
and dust grain opacities (e.g.. Pollack et al. 1994). The code has been tested on a variety of 
test problems (Boss & MyhiU 1992; MyhiU & Boss 1993). 

The Poisson equation for the cloud's gravitational potential is solved by a spherical 
harmonic expansion {Yim) including terms up to Nim = 32. The computational grid consists 
of a spherical coordinate grid with Nr = 200, Ng = 22 for 7r/2 > ^ > (symmetry through 
the midplane is assumed for tt > ^ > 7r/2), and A^^ = 256 for 27r > > 0, i.e., with no 
symmetry assumed in (p. The radial grid contracts to follow the collapsing inner regions and 
to provide sufficient spatial resolution to ensure satisfaction of the four Jeans conditions for 
a spherical coordinate grid (Truelove et al. 1997; Boss et al. 2000). The innermost 50 radial 
grid points are kept uniformly spaced during grid contraction, while the outermost 150 are 
non-uniformly spaced, in order to provide an inner region with uniform spatial resolution 
in the radial coordinate. The (j) gi"id is uniformly spaced, whereas the 9 grid is compressed 
toward the midplane, where the minimum grid spacing is 0.3 degrees. 



Tables 1 and 2 list the initial conditions for the models. The initial models have Gaussian 
radial density profiles (Boss 1997) of the form 



where po = 2.0 x 10~^^ g cm~^ is the initial central density. The prolate clouds with central 
densities 20 times higher than a reference boundary density have ra — 1.16R and Vb — rc — 
0.580i?, where R is the cloud radius, yielding a axis ratio of 2 to 1. The oblate clouds with 
the same degree of central concentration have Va = ri, = 1.16i? and Tc = 0.580i?. For 100 to 
1 density contrasts, the prolate clouds have ra — 0.932i? and Vb — Tc — 0.466i?, while the 
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Initial Conditions 




(1) 
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oblate clouds have Va — rb — 0.932R and Tc = 0A66R. Random numbers {ran{x, y, z)) in the 
range [0,1] are used to add noise to these initial density distributions by multiplying pi from 
the above equation by the factor [1 + 0.1 ran{x,y, z)]. The cloud radius is R — 1.0 x 10^^ 
cm w 0.032 pc for all models. 

The cloud masses are 1.5 Mq and 2.1 M©, respectively, for the prolate and oblate clouds 
with a 20:1 density ratio, and 0.96 Mq and 1.5 Mq, for the prolate and oblate clouds with 
a 100:1 density ratio. With an initial temperature of 10 K, the initial ratio of thermal to 
gravitational energy is ai = 0.39 for the prolate clouds and 0.30 for the oblate clouds with 
the 20:1 initial density ratio, while ctj = 0.55 for the prolate clouds and 0.39 for the oblate 
clouds with the 100:1 initial density ratio. 

Solid body rotation is assumed, with the angular velocity about the z axis (short axis) 
taken to be Qi = 10^^^, 3.2 x 10~^^, or 10~^^ rad s^^. These choices of Qi result in initial ratios 
of rotational to gravitational energy varying from (3i = 1.2 x 10^^ to 0.015 for the prolate 
clouds and Pi = 1.1 x 10~^ to 0.013 for the oblate clouds. These choices are all consistent 
with observational constraints on the densities, shapes, and rotation rates of dense cloud 
cores (e.g., Myers et al. 1991; Goodman et al. 1993; Ward-Thompson, Motte, & Andre 
1999; Jones, Basu, & Dubinski 2001; Curry and Stabler 2001). 

As in the previous three-dimensional models (Boss 1997, 1999, 2002, 2005, 2007), the 
effects of magnetic fields are approximated through the use of several simplifying approxima- 
tions regarding magnetic pressure, tension, braking, and ambipolar diffusion (see Boss 2007 
for a derivation of these approximations). All models assumed an ambipolar diffusion time 

scale tad = 10 tff, where the free fall time tff — (37r/32G'po)^^^ = 4.7 x 10^ yr for a central 
density po = 2.0 x 10^"*^^ g cm^^. The reference magnetic field strength is assumed to be 
-Boj = 200/iG. The magnetic braking factor f^b is taken to be either 0.0001 or 0.001. Based 
on the models of Basu & Mouschovias (1994), Boss (2007) estimated that f^b ~ 0.0001. 
The models with fmb = 0.001 are thus intended to attempt to overestimate the effects of 
magnetic braking. 

Note, however, that the calculations of Basu & Mouschovias (1994) stopped at central 
densities of 3 x 10^ cm~^, i.e., before the clouds became optically thick, whereas the present 
models are continued well in the optically thick regime. Hence the implicit assumption that 
the trends found in the Basu & Mouschovias (1994) models and used by Boss (2007) to 
derive the fmb approximation should continue indefinitely may not be warranted, though 
the trends do persist over the previous six orders of magnitude increase in central density 
of their models (see Figure 7 of Basu & Mouschovias 1994). The magnetic braking studied 
by Basu & Mouschovias (1994) is similar to the magnetic braking produced by disk jets and 
outfiows in the models by Banerjee & Pudritz (2006), Hennebelle & Fromang (2008), and 
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Hennebelle & Teyssier (2008). A superior treatment of magnetic braking beyond the fmb 
approximation of Boss (2007) will require a true MHD code. 

With the field strength B^i — 200//G, the prolate or oblate clouds with 20:1 density con- 
trasts have initial ratios of magnetic to gravitational energy of 7j = 0.58 or 0.43, respectively. 
For density ratios of 100:1, 7^ = 0.81 for the prolate clouds and 0.57 for the oblate clouds. 
The mass to flux ratio of these clouds is less than the critical mass to flux ratio, making all 
the clouds formally magnetostatically stable and hence magnetically subcritical. Protostel- 
lar collapse cannot occur until ambipolar diffusion leads to sufficient loss of magnetic field 
support for sustained contraction to begin. 

4. Results 

Tables 1 and 2 list the initial conditions as well as the basic outcome of each model, 
namely the final time to which the cloud was advanced tf (in units of the initial cloud free 
fall time) and whether the cloud collapsed to form a quadruple system (Q), binary (B), 
binary-bar (BB), single-bar (SB), ring (R), or did not collapse (NC). For convenience, the 
results of the corresponding models by Boss (2002) without magnetic braking are shown as 
well. 

Figure 1 shows the initial equatorial density distribution for the prolate cores with a 20:1 
density contrast. Given the stability of the initial models, evolutions consist of the clouds 
oscillating about the initial configurations, waiting for sufficient time to elapse for ambipolar 
diffusion to reduce the magnetic field support enough to allow collapse to proceed, as in 
the previous magnetic cloud models (e.g.. Boss 1997). Figure 2 shows the result for model 
P2BB, which collapsed to form a quadruple protostellar system, though in this case with an 
additional central density maximum. Without magnetic braking, this core collapsed to form 
a binary protostar system, so in this case magnetic braking has led to an increased degree of 
fragmentation. The off-axis clump evident in Figure 2 at 9 o'clock has a maximum density of 
2.5 X 10^^^ g cm~^ and an average temperature of 20 K. Considering regions with a density at 
least 0.1 that of this maximum density, the clump's mass is 2.6 Jupiter masses, greater than 
the Jeans mass at that density and temperature of 1.9 Jupiter masses. The ratio of thermal 
to gravitational energy for the clump is 0.49, showing that it is gravitationally bound. The 
ratio of rotational to gravitational energy is 0.39, so the clump is in rapid rotation. These 
fragment properties are quite similar to those found in the Boss (2002) models. 

Model P2BA behaved in much the same way as P2BB, even with stronger magnetic 
braking (/^^ = 0.001 for P2BA compared to f^b = 0.0001 for P2BB). 
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Figure 3 shows the outcome of model P2BD, which started with a lower initial angular 
velocity than model P2BB but was otherwise identical. In this case, a well-defined binary 
protostar system forms. A similar outcome resulted for model P2BC with a higher degree 
of magnetic braking. With even lower initial angular velocity than P2BC and P2BD, model 
P2BE collapsed to form a binary-bar system (Figure 4), i.e., a binary with its members 
connected by a bar of rotating gas. In the case of model P2BE, there are two density 
maxima in the bar near the center of the system, possible evidence for further fragments, 
though none of the density maxima are as well-defined as those of the binary in Figure 3. In 
contrast, when magnetic braking was neglected (Boss 2002), a core identical to model P2BE 
fragmented into a well-defined binary system, so in this case magnetic braking has reduced 
the degree of fragmentation. 

The prolate core models with 100:1 initial density contrasts all behaved in the same 
manner and formed single central protostars embedded in bar-like structures. Figure 5 shows 
the result for model PIBB, termed a single-bar. Figure 6 shows the equatorial temperature 
contours for this model, emphasizing the highly non-uniform temperature distribution that 
results from including radiative transfer effects (e.g.. Boss et al. 2000). The temperature 
field is strongest in regions where the infalling gas is forming shock-like density corrugations, 
yielding an x-shaped pattern in the cloud's midplane. The formation of single-bars in all five 
of the prolate 100:1 models compared to the formation of binaries in the corresponding Boss 
(2002) models shows that when the initial cores are highly centrally condensed, magnetic 
braking is able to frustrate fragmentation. 

Figure 7 shows the initial conditions for the oblate cores with 20:1 initial density con- 
trasts, while Figure 8 depicts the outcome of model 02BB: a well-defined ring. While the 
ring shows no particular tendency to fragment over the time scale of the calculation, such 
a configuration is expected to fragment eventually. With the exception of model 02BA, 
which did not collapse significantly, all of the oblate cores collapsed to form rings, for both 
20:1 and 100:1 density contrasts, though the rings were more pronounced for the cores with 
higher initial rotation rates, such as 02BB in Figure 8. In comparison, the corresponding 
Boss (2002) models formed a combination of rings or rings which showed a tendency to 
fragment into quadruple systems. Hence the oblate cloud models indicate that the inclusion 
of magnetic braking had only a mild tendency to inhibit their fragmentation. 

Tables 1 and 2 show that the prolate clouds tended to take considerably longer to 
undergo collapse than the oblate clouds, implying a considerably longer period in the pre- 
collapse, quasi-equilibrium phase where thermal and magnetic support dominate. As a 
result, the effect of magnetic braking through the fmb approximation should be stronger in 
the models which took the longest time to collapse, i.e., the models in Table 1 with Po/pr — 
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100 : 1 and tf/tffr^9. Prolonged magnetic braking will tend to suppress rotationally-driven 
fragmentation, consistent with the formation of single-bars for the prolate models in Table 1 
with tf/tff ~ 9. Model P2bb, for example, lost about 0.4% of its total angular momentum 
during its evolution to tf/tff = 4.683, compared to a loss of about 1% for model PlBb within 
tf/tff = 9.071. Note though that for the oblate models in Table 2, no such effect is evident, 
perhaps because all of those models collapsed within tf/tff < 5. In fact, a comparision of 
the corresponding oblate models shows that model 02bb lost about 0.4% of its total angular 
momentum during its evolution to tf/tff = 2.456, compared to a loss of about 1% for model 
OlBb within tf/tff — 4.636. Both models 02BB and OIBB collapse to form rings. Their 
percentage angular momentum losses are identical to those for models P2BB and PIBB, 
which formed a quadruple and single-bar, respectively, implying that the initial cloud shape 
and degree of central concentration do have an important effect on the fragmentation process 
of magnetic clouds. 



5. Conclusions 

These pseudo-magnetohydrodynamics calculations have explored the possibly deleteri- 
ous effects of magnetic braking on protostellar fragmentation, an issue explored by Hosking 
& Whitworth (2004). A degree of inhibition of fragmentation caused by magnetic braking of 
both prolate and oblate, dense cloud cores has been identified in these models through a di- 
rect comparison with otherwise identical models calculated by Boss (2002) without magnetic 
braking. Nevertheless, a cursory examination of the figures and tables reveals that there is a 
still a large portion of initial conditions space that appears to be permissive of fragmentation 
of magnetic clouds when the approximate effects of magnetic pressure, tension, and braking 
are all included. The present models thus suggest that binary and multiple stars may well 
form from the collapse and fragmentation of magnetic, as well as non-magnetic, dense cloud 
cores, though perhaps not quite so readily. 

Given the approximate nature of the present calculations (Boss 2007), it will be impor- 
tant to try to confirm these results with a true magnetohydrodynamics code. This could be 
accomplished by adding a numerical solution of the magnetic induction equation to the Boss 
& Myhill (1992) code. Alternatively, these calculations could be repeated using publically 
available MHD codes, such as the FLASH adaptive mesh refinement code (e.g., Duffin & 
Pudritz 2008), though FLASH docs not at present include Eddington approximation radia- 
tive transfer, unlike the Boss & Myhill (1992) code. Attempting such true MHD calculations 
stands as a challenge for future work on the question of protostellar collapse and fragmen- 
tation. 
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Table 1. Initial conditions and results for initially prolate cores. In this table and the 
following, po/ pR is the initial ratio of the central to the boundary density, the units for the 
initial angular velocity fij are rad s^^, and the magnetic braking factor f^b is 
dimensionless. The final times tf are given in units of the initial free fall time 
tff = {3n/32Gpoy/^ = 4.7 x 10^ yr. Q denotes a core that collapses to form a quadruple 
system, B denotes a binary outcome, BB a binary-bar, and SB a single-bar. The results 
obtained by Boss (2002) for identical cores but without magnetic braking are shown as 

well. 



model 


Pol PR 




fmb 


tf/tff 


result 


Boss(2002) 


P2BA 


20:1 


1.0 X 10-13 


0.001 


4.649 


Q 


B 


P2BB 


20:1 


1.0 X 10-13 


0.0001 


4.683 


Q 


B 


P2BC 


20:1 


3.2 X 10-1^ 


0.001 


4.602 


B 


B 


P2BD 


20:1 


3.2 X 10-1^ 


0.0001 


4.603 


B 


B 


P2BE 


20:1 


1.0 X 10-1^ 


0.001 


4.833 


BB 


B 


PIBA 


100:1 


1.0 X 10-13 


0.001 


9.065 


SB 


B 


PIBB 


100:1 


1.0 X 10-13 


0.0001 


9.071 


SB 


B 


PIBC 


100:1 


3.2 X 10-1^ 


0.001 


8.942 


SB 


B 


PIBD 


100:1 


3.2 X 10-1"! 


0.0001 


8.945 


SB 


B 


PIBE 


100:1 


1.0 X 10-1^ 


0.001 


8.918 


SB 


B 
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Table 2. Initial conditions and results for initially oblate cores. R denotes a ring formed, 
while NC means no significant collapse occTirred. 



model 


Pol PR 




fmb 




result 


Boss(2002) 


02BA 


20:1 


1.0 X 10-^3 


0.001 


4.565 


NC 


R 


02BB 


20:1 


1.0 X 10-^3 


0.0001 


2.456 


R 


R 


02BC 


20:1 


3.2 X 10"^^ 


0.001 


2.127 


R 


Q 


02BD 


20:1 


3.2 X IQ-^* 


0.0001 


2.134 


R 


Q 


02BE 


20:1 


1.0 X 10-14 


0.001 


2.109 


R 


Q 


OIBA 


100:1 


1.0 X 10-13 


0.001 


4.656 


R 


R 


OIBB 


100:1 


1.0 X 10-13 


0.0001 


4.636 


R 


R 


OIBC 


100:1 


3.2 X 10-14 


0.001 


4.521 


R 


R 


OIBD 


100:1 


3.2 X 10-14 


0.0001 


4.495 


R 


R 


OIBE 


100:1 


1.0 X 10-14 


0.001 


4.422 


R 


Q 
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RHOMAX= -17.7 CONDIF= 0.1 R= O.lOE+18 




Fig. 1. — Initial equatorial density contours for prolate core models with 20:1 density con- 
trasts. Maximum density is 2.0 x 10~^^ g cm~^. Contours represent changes by a factor of 
1.3 in density. Region shown is 1.0 x 10^^ cm in radius. 
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Fig. 2. — Equatorial density contours for model P2BB at a time of 4.683 tjf. Maximum 
density isl.OxlO~^^g cm~^. Contours represent changes by a factor of 1.3 in density. Region 
shown is 4.0 x 10^^ cm in radius. A possible quintuple protostellar system has formed: a 
central density maximum surrounded by four density maxima. 
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Fig. 3. — Equatorial density contours for model P2BD at a time of 4.603 tff. Maximum 
density is 1.3 x 10~^^ g cm~^. Contours represent changes by a factor of 1.3 in density. 
Region shown is 4.0 x 10^^ cm in radius. A binary protostellar system has formed. 
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Fig. 4. — Equatorial density contours for model P2BE at a time of 4.833 tff. Maximum 
density is 1.6 x 10"^^ g cm"^. Contours represent changes by a factor of 1.3 in density. 
Region shown is 4.0 x 10^"^ cm in radius. A binary-bar system with four different local 
density maxima has formed. 
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Fig. 5. — Equatorial density contours for model PIBB at a time of 9.069 tff. Maximum 
density is 6.3 x 10~^^ g cm"^. Contours represent changes by a factor of 1.3 in density. 
Region shown is 4.0 x 10^^ cm in radius. A single-bar system has formed. 




Fig. 6. — Equatorial temperature contours for model PIBB at a time of 9.069 tff. Maximum 
temperature is 25 K. Contours represent changes by a factor of 1.3 in temperature. Region 
shown is 4.0 x 10^^ cm in radius. 
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RHOMAX= -17.7 CONDIF= 0.1 R= O.lOE+18 




Fig. 7. — Initial equatorial density contours for oblate core models with 20:1 density con- 
trasts. Maximum density is 2.0 x 10~^^ g cm~^. Contours represent changes by a factor of 
1.3 in density. Region shown is 1.0 x 10^^ cm in radius. 
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Fig. 8. — Equatorial density contours for model 02BB at a time of 2.456 tff. Maximum 
density is 2.0 x 10~^^ g cm~^. Contours represent changes by a factor of 1.3 in density. 
Region shown is 2.3 x 10^'' cm in radius. A ring has formed: a strong density minimum 
occurs at the center of the core. 



